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^H Abstract 

O 

^NJ The recent measurements on Rk and Rt^ imply that there exists a possible violation of the 

Q^ leptonic flavor universality which is one of the cornerstones of the standard model. It is suggested 



< 



that a mixing between sterile and active neutrinos might induce such a violation. In this work 

^_^ we consider the scenarios with one or two sterile neutrinos to explicitly realize the data while the 

Q^ constraints from the available experiments have been taken into account. Moreover, as indicated in 

^^ literature, the deviation of the real PMNS matrix from the symmetric patterns may be due to a fJ,—T 

II asymmetry, therefore the measurements on RD{Ds)eii = r(D(-Ds) — )■ e'^Ve)/'^{D{Ds) — ^ f-^i^ii) and 

1— i Rd(Ds)ut — ^{D{Ds) -^ f'^'^n)/'^{D{Ds) -^ f-^Tr) (and for some other heavy mesons B^ and Be 

^^ etc.) may shed more light on physics responsible for the violation of the leptonic flavor universality. 

^_J The data of BES III are available to test the universality and that of the future charm-tau factory 

■^ will provide more accurate information towards the aspect, in this work, we will discuss RD{Ds)ea 

^^ and RD(Ds)fiT ill a-H details and also briefly consider the cases for B and Be- 
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PACS: 14.60.Pq Neutrino mass and mixing 



I. INTRODUCTION 

The property that couphngs of leptons to gauge bosons are independent of lepton flavors 
is named as lepton flavor universality (LFU), which is embedded in the Standard Model 
(SM) and any violation of LFU may be induced by new physics beyond the SM. Taking 
leptonic decays of W boson W -^ Ivi (/ = e, /i, r) as an instance, ratios of the branching 
fractions measured at LEP 11 [Ij are 

B{W -^ iiV^,)/B{W -^ eue) = 0.997 ± 0.021, (1) 

B{W -^ TVr)/B{W -^ eVe) = 1.058 ± 0.029, (2) 

B{W -^ TVr)/B{W -^ iw^) = 1.061 ± 0.028. (3) 

which demonstrate the lepton universality at the level of 2.9%. Besides, leptonic decays of 
mesons also provide a possibility to test LFU where the uncertainties in the hadronic sector 
are cancelled. In literature the ratios are suggested to be measured 

where P = n, K, D, Dg, B, B^. and a, /3 = e, /i, r. In order to clearly demonstrate deviation 
of the measured value R^p^a from the SM predictions -Rp^^, a parameter Arp^^ is defined 
as 

Arp„;3 = -^ - 1. (5) 

To accommodate non-zero Arpap, one may invoke two different mechanisms [2]: 

• Introducing new Lorentz structure in the four-fermion interaction; 

• Modifying the Wli^i vertex by corrections to lepton mixing. 

For the first category, SM may be extended to new physics beyond standard model (BSM) 

which includes charged Higgs, for example in the supersymmetry (SYSY) [3] or two-Higgs- 
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doublet models [1]. In the SM, the violation of LFU are estimated as |Arp^^| = 0{-^ 
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which is too small to be accounted, while for the first category, the correction to the Wav^ 
vertex emerges at loop level where new physics BSM particles exist is of order 0(^ x ^^), 
which is greatly suppressed by the new physics scale Anp [Sj. 



For the second category, the couphng Wava is modified by breaking the unitarity of the 
3x3 lepton mixing matrix. In the type-I seesaw [B], sterile neutrinos are introduced and 
the 3x3 mixing matrix could be non-unitary to the level of 0{M^/M^) [7] [5], where ttid 
and mji are the Dirac and Majorana mass matrices, respectively. 

Introducing sterile neutrinos can provide natural interpretations of some anomalies ob- 
served in recent experiments. For instance, the LSND p], the MiniBooNE [10], the reactor 
[TT] anomalies and as well as the gallium anomaly of the GALLEX [12] and SAGE [13] ex- 
periments. Furthermore, a careful analysis of the cosmological data show that the effective 
number of neutrino species is larger than 3 [14J , which might also hint existence of additional 
one or two neutrino species besides the three active ones. Therefore a study on violation of 
LFU would be interesting because it may reveal possible new physics mechanisms beyond 
the SM at the lepton sector. 

As indicated in literature, the deviation of the real PMNS matrix from the symmetric 
patterns may be due to a /i — r asymmetry. Since the kinematic constraints, K and vr cannot 
decay into rz/, the measurements on RD{Ds)efi = r(-D(-Ds) — t- e~^h'e)/^{D{Ds) — )• yU+z/^) and 
Rd{Ds)plt = r(-D(-Ds) — ^ /^^^M)/r(-D(-Ds) — !■ fJ'~^Tr) (and some other heavy mesons S^ and Be 
etc.) may shed more light on the physics responsible for the violation of the leptonic flavor 
universality. The data of BES III are available to test the universality and measurements 
at the future charm-tau factory will provide more accurate information towards the aspect. 
Indeed, the experimental results of NA62 Collaboration [T3] determine Arxefi = (4±4) x 10~^ 
and Ar^e/. = (-4 ± 3) x 10^^ ^ whereas the BES data tell us Aro^^r = 0.256+°;^?° which 
is rather large. 

The presence of sterile neutrinos changes the leptonic decay widths of pseudoscalar mesons 
and then Arp^^ may be enhanced or suppressed. In this work we start with this motivation 
and derive the analytical formulas of Apq,/3 for 3+1 and 3+2 scenarios in the next section. 
Especially in this work, we will discuss RD(Ds)efi and -Rd(Ds)/xt in all details and also briefly 
consider the cases for B^ and Be- Numerical analyses are made and possible experimental 
measurements on Arpa/3 at BES and future charm-tau factory are discussed. Then we 
present some discussions in the last section. 



II. Rpap WITHIN SM 

In the SM, the leptonic decay widths of pseudo-scalar mesons P — )■ av^ (a = e, /i or r 
and P = IT, K, D, Dg, B or Bf) are given by 

r|".... = g|v^.'lVI-P<(i-^)', (6) 

where Gp is the Fermi constant, ma is the lepton mass, Vqqi the element of Cabibbo- 
Kobayashi-Maskawa (CKM) mixing matrix corresponding to the constituents in the meson 
P whose mass is irtp and decay constant fp and neutrino masses are neglected. As well 
known, these processes are helicity suppressed. 

With the definition of Fp^Q,^^ in Eq. (p|, one has 

^''-^-[mj [ml -ml) ' ^^^ 

If the QED corrections are considered, there will be a factor (1 + Sqed) multiplying to the 
left-handed side of Eq. ([T]). An approximate estimate of electrodynamical correction to 
P — !■ aua is about 1I7 x 2^ ~ 10~'^ suppression. For the leptonic kaon decay, due to the 
inner bremsstrahlung Ki2^ process which is included by definition into Rk, the calculation 
indicates 6qed = (-3.78 ± 0.04)% [IB]. 

As suggested, existence of sterile neutrinos changes the picture, so in the next section, 
we will study how they contribute to P — )■ oz/q,. 

III. LEPTON MIXING WITH PRESENCE OF STERILE NEUTRINOS 

Now let us give a general description of the involvement of sterile neutrinos and take 3-1-1 
as an illustration. In the SM, the sterile neutrinos do not directly couple to PF-boson, thus 
without a mixing between sterile neutrinos with active ones, the weak interaction can be 

written as 



(e,/i,r)7'^(l-75)(?7pMNS,0) 



1^2 



W,, (8) 



where Piii = 1,2, 3) are the active neutrinos and z/^ is a sterile neutrino (or several), t/pMNS is 
the regular PMNS lepton mixing matrix and here denotes a 3x1 matrix, thus (f/pMNS)0) 



is a 3 X 4 matrix and in this scenario, unitarity is no longer existing. When the mixing 
between z/^ and z/j is introduced, the matrix becomes (f/pMNS " cose, sine) and the couphng 
vertex is 



(e, ju, r)7^(l - 75) (f/pMNS ■ cos e, sin e) 



T^2 



W, 



M' 



(9) 



then the sterile neutrino participates in the weak interaction^. Meanwhile this mixing also 
induces a flavor-changing neutral current (FCNC) as 



(z/i,z/2,i/3,i^.)7^(l-75)(l,0) 



T^2 



Zf, -^ (z/1, z/2, z/3, z^si)7^(l - 75)^' 



1^2 



^A" 



(10) 



where V^' is a 4 x 4 matrix and resulted in from a matrix product of (4 x 3) ® (3 x 4) i.e. 
(f^PMNS ■ cos e, sin e)''"(f/pMNS ■ cos e, sin e). This non-unitary mixing matrix and sizable sterile 
neutrino masses can also change the generation number which is accurately measured by 
the LEP experiment. Therefore the mixing is rigorously constrained by the LEP data. 

Once sterile neutrinos are introduced, they do mix with active ones, the decay width 
obtained above (|6l) will be modified into 
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where tti^j is the sterile neutrino masses and A^ is the total number of neutrino mass eigen- 
states, for instance, assuming two sterile neutrinos, i.e., the so-called 3 + 2 scenario [T7] . 
A^ = 3 + 2 = 5. The A(mp, m^, m^J is defined as 



A(mp, m^, m^j) = \jvn\, + m^ + m).^ — 2vr?pV[i\ — 2m 
Now we can obtain the general expression of Apq,^ 



pm% - mlmli. 



Arp. 
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1, 
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^ After the replacement, 1^1,1^2,1^3, ^s are slightly changed and constitute the real mass eigenstate. 



where R^a is obtained in Eq. (P7| and 
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y + E,=i l^/3j+3p(rw^ + m2^.) (m| - m2 - m^^.) A( 
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(14) 



A. 3+1 scenario 

In order to get the 3x4 mixing matrix describing the sterile-active neutrino mixing 
in charged- current and neutral-current interactions, an easy means is to construct a 4 x 4 
matrix in a regular way and then remove away the last row of the 4x4 matrix to get the 
required 3x4 matrix. The 3x5 matrix for the 3+2 scenario will be obtained in a similar 
way. 

First we deal with the case by adding one sterile neutrino into the game, i.e., the 3 + 1 
scenario, in which the lepton mixing matrix can be parameterized with 6 mixing angles and 
3 phases in a form 



U: 



3+1 



R. 



34l,t734, 



'^34)-R24(6'24, 0)i?i4(6'i4, Sm)R23{023, 0)-Ri3(6'i3, '^13)-Ri2(6'i2, 0), 



(15) 



where Rij is a 4 x 4 matrix describing the rotation in the i — j plane and here presents i?34 
as an instance, 



-^34(^34,(^34) 



/l 

1 







COS 6^34 



\0 -sin ^346 



1034 






sin 6'34e" 
cos 6* 



\ 



1634 



(16) 
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It is obvious that the multiplication R23R13R12 = U^ is the usual Pontecorvo |18j-Maki- 
Nakawaga-Sakata [19j (PMNS) matrix modified by adding a trivial fourth column and row. 
In order not to be contradict with the experimental measurements of available three- 
neutrino oscillations and the data of LEP, the mixing between sterile and active neutrinos 
must be very small. And here an assumption is made that the added sterile neutrino does 
not distinguish between the three active ones, i.e., the mixing angles 6*34 = 6^24 = 6*14 = ei. 



In Eq. (14) only squared mixing elements \Uai\ appear in Arpai3 while neglecting higher 



order powers of ei, one obtains 

\Uai\^ = |?7°ipcos^ei,(a = e,fi,T;i = 1,2,3); 

17-7-12 -2 lr7-|2 -2 2 \TT \2 -2 4 /iTN 

\Ue4\ = sm ei; |ty^4| = sm ei cos ei;\Ur4\ = sm ei cos ei, (17j 

where we ignore the subscript "3+1" in t/3+1 and the f/^j refers to the elements of the unitary 
f/o (3+0) with J2Li K? = 1 for a = e,/i,r. 

Then with these mixing matrix elements, one can obtain R^p^a in the 3+1 scenario 

cos^eim^(m^ - mlY + \Ua4\'^iml + mfj (mp - m^ - m^^) A(mp, m^,m^J 



i?l 



COS' 



eim^(mp — m^)^ + |t/g4p(m^ + m^J (rrip — "^^ — ^^1) ^{^'p-, ^% ^li) 



(18) 



which is substituted into Eq. (13) to obtain Arpa/B under the 3+1 scenario. 



B. 3+2 scenario 

Turning to the 3+2 scenario in which two sterile neutrinos are introduced, the procedure 
is similar to that for the 3+1 case. The mixing angles between the second sterile neutrino 
and active neutrinos are denoted as €2- Then we summarize the results below: 

\Uai\'^ — \U^i\'^cos'^ eicos^ 62, {a = e,fi,T;i = 1,2,3); 
\Ue4\ = sin ei cos €2; jt/esl = sin 62; 

I7-7- |2 -2 2 2 17-7- |2 • 2 2 

I'-^tj.i] = sm ei COS ei cos 62; Itv^sl = sm e2Cos 62] 

\Ur4:\ ~ sin ei cos ei cos €2; |f/r5| = sin €2 cos 62- (19) 

With these elements one can obtain Arpap with the 3+2 scenario. 

C. Constraints set by the FCNC 

With the scenarios discussed above, the neutrino neutral current (NC) interaction will 
be modified. As discussed above, the number of neutrinos is 2.984 ± 0.008 by fits to LEP 
data. Assuming the small deviation from 3 to be caused by mixing between active and 
sterile neutrinos, this value will cast rigorous constraints on the mixing parameter ei_2. An 
estimate is: ei < 0(5 x 10"^) for the 3+1 scenario and e\ + e\< O(10~^) for the 3+2 case. 



Scenario 


\Uf.e\ 


\Ure\ 


I^tmI 


3+1 


ei < 1.75 X 10-3 


ei < 3.70 X 10-2 


ei < 3.599 X 10"^ 


3+2 


ef + ei < 3.05 x 10^^ 


el+el< 1.37 X 10-^ 


ef + ei < 1.295 X 10-^ 



TABLE I: Constraints to new added mixing parameters ei^2 from the experimental limits of LFV 
processes. 



A preliminary result about sum X]i=i Wei\ = 0.994 ± 0.005 at 90% confidence level [20] , 
which signifies the non-unitarity of the 3x3 active neutrino mixing matrix to be at level 
< O{10~^) [21] • The non- vanishing terms Uai3 = J2i=4 ^ai^iP ('^ 7^ P) result in the tree level 
FCNC interaction, which can induce the low-energy lepton flavor violation (LFV) processes. 
And these LFV precesses are proportional to the value of the non-vanishing X]i=4 ^ai^iis- 
The experimental bounds on these LFV interactions can be transformed to constraints on 
the new mixing angles ei^2- For instance, constraints to |f/a^| from several LFV processes 
Ref. [221123]: 

\U^e\ < 3.05 X 10-^ {B{fi- -^ e-e+e-) < 1.0 x 10-^^) , 
\Ure\ < 1.37 X 10"^ {B{t- -^ e-e+e-) < 3.6 x 10"^) , 
\Ur^,\ < 1.295 X 10-3 {B{r- -^ fi-fx+f^-) < 3.2 x 10-^) . 

Considering these bounds and the elements derived above, with reasonable approximations 
we can obtain constraints on ei_2 which are shown in Table [T] 

IV. NUMERICAL ANALYSES 



In this section we numerically evaluate the total contributions from SM and new physics 
BSM to ArpaiB in the context of 3+1 and 3+2 scenarios to account for the observed data. The 
analyses are about the leptonic decays of pseudo-scalar mesons P = tt,K, D, Dg, B, Be- The 
evaluated branching ratios of the leptonic decays are listed in Table [IT] and the corresponding 
Rpap and Arpajs are shown in Table III 



p 


mass (MeV) 


B{P -^ ei^e) 


BiP ^ fii^^) 


B{P -^ TVr) 


-K 


139.57018 ± 0.00035 


(1.230 ±0.004) X 10-^ 


(99.98770 ± 0.00004)% 


- 


K 


(493.677 ±0.016) 


(1.581 ±0.008) X 10-5 


(63.55 ±0.11)% 


- 


D 


(1869.62 ±0.15) 


< 8.8 X lO-'^ 


(3.82 ±0.33) X 10-4 
[(4.15l[J:22) X 10-4] 


< 1.2 X 10-3 
[(1.10 ±0.06) X 10-3] 


Ds 


1968.49 ±0.32 


< 1.2 X 10"^ 


(5.90 ±0.33) X 10-3 
[(5.50lO:i) X 10-3] 


(5.43 ±0.31)% 
[(5.36l°j4) X 10-2] 


B 


5279.25 ±0.17 


< 9.8 X 10-' 

[[(1.1 ±0.2) X 10-11]] 


< 1.0 X 10-6 

[[(4.5 ±1.0) X 10-^]] 


(1.65 ±0.34) X 10-4 

[(0.796;n87) X 10-4] 
[[(1.0 ±0.2) X 10-4]] 


Be 


6277 ± 6 


- 


- 


- 



TABLE II: The experimental values or bounds of masses and branching ratios for P = 
vr, K, D, Ds, B, Be taken from PDG [22j. The corresponding theoretical predictions are in (double) 
square brackets which are from Ref. ([ 



p 


ocxp 


«, 


^rpef. 


TT 


(1.230 ±0.004) X 10-4 


1.234 X 10-4 


(-3.241 ±3.241) X 10-3 


K 


(2.488 ±0.013) X 10-5 


(2.472 ±0.001) X 10-5 


(6.472+5-668) ^ 10-3 


P 


RT,r 


R'pllr 


Arp^^ 


D 


> 0.291 


n Q 77+0. 043 
^•"-" '-0.038 


> -0.308 


Ds 


n inq+o-oi3 
u.iuy_o,oi2 


(8.65+J:^|) X 10-2 


0.25610;430 


B 


< 7.6 X 10-3 


(4.5+2-4) X 10-3 


< 1.6 


Be 


- 


(4.18+n^) X 10-3 


- 



TABLE III: The current experimental measurements and SM prediction of Rpefj, [25j, Rpfj^r and 
the corresponding Arpai3- The SM predictions include the uncertainties from electromagnetic cor- 
rections and as well as the uncertainties due to CKM mixing matrix elements and decay constants. 
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FIG. 1: The allowed parameter space of ei and m^i for Arj^eu and Arx, 
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A. 3+1 scenario 

In this subsection, we make a numerical analysis on the breaking of the lepton universality 
in different decay processes, within the 3 + 1 scenario. Ar^g^ and Arxefi in the parameter 
space of mixing parameter ei and sterile neutrino mass nisi are shown in Fig. [TJ The central 
value of Arpe/x is denoted by the solid line while its 1 — cr region is enclosed by the dashed 
lines. In the left panel of Fig. [Tlfor Ar^g^, to accommodate the experimental results, the 
mass of the sterile neutrino has to be larger than 250 MeV, thus the final state with Ug is 
kinematically forbidden. For ArKen, the allowed parameter space of ei and rrisi is larger, 
which covers ~ 400 MeV and ~ 0.35, respectively. Also, the region of rusi > ttlk, i^s 
does not show up in the final state and the region with relatively larger ei needs to be ruled 
out for its failure to be reconciled with the LEP data constraint. 

Besides, as shown in Fig. [2| the 3+1 scenario fails to provide a common parameter space 
of ei and rrisi to saturate the experimentally measured Ar^e^ and Ar^e^ simultaneously, 
namely there does not exist a solution within la tolerance. 

Since both Ds to fiu and tu have been experimentally measured, thus Ar^)^^,- is obtained. 
Lack of experimental data on the decay rates of leptonic decays oi D, we cannot determine 
Ato^t yet, for an illustration, we set the ArD^^T as 10^"*^, 10^^, 10~'^, 10""^, 10~^ to get a sense 
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FIG. 2: Comparison of ArT^e/i and Arxefi in the parameter spaces of mixing angles ei and the 
sterile neutrino mass nisi- 

about the dependence of the Ars on ei and rrisi- The results are shown in Fig. |3} 

Form the figures of Ar^e^, Ar^i^e/x and Ar^^T- in Fig. [3l it is obvious that non-zero Ar 

demands non- vanishing ei. Moreover, it is noted that for fixed rrisi (ei), the smaller Ar is, 

the smaller ei {nisi) would be. 

In (d) of Fig. pi the lower bound of Ar^^^r = 0.25Qt.o'^i^ does not appear, as our 

analyses indicate that Ar^^nr cannot be negative in the 3+1 scenario. From this diagram, 

it is obvious that within 1 — a range of Ar^)^^^, the particular values ei = and (rrisi = 0) 

are not excluded and the vanishing sterile-active neutrino mixing signifies that the lepton 

fiavor universality holds. Thus to make a decisive judgement more accurate measurements 

are needed. 

As discussing leptonic decays of B and Be mesons, because of lack of experimental data, 

we take several values for Ar to illustrate its dependence on the parameters as shown in 

Fig. i 
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FIG. 3: Aroeti, ^rD.e/i, ^rofir, ^ro^^ir vs ei and nisi. 
B. 3+2 scenario 

In this subsection, we numerically analyze the lepton universality with two sterile neu- 
trinos, i.e., the 3+2 scenario. 

Even though the errors are still large, Ar^(x)e/i have been set, thus we first present Arjre/i 
and ArKe/i in the same graph Fig. [Sl where the horizontal solid line corresponds to the 
central value and dashed lines enclose the 1 — cr range of ArxetJ. whereas the perpendicular 
ones are for Ar^e^t- 

The cross region satisfies both Ar^re/x and Ar^-g^. In our analyses, we let the mixing angles 
ei and €2 vary within (0, 3 x 10"^) and (0, 5 x 10~^), respectively, while the sterile neutrino 
masses nisi ^ (0, 140) MeV and ms2 G (0, 500) MeV. It is noted that for such parameter 
ranges, there exist solutions to accommodate both Ar^e^ and Arxeti- Concretely, the red 
points which correspond to the values calculated within the parameter ranges fall in the 
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FIG. 4: Arsefj,, Ar^^e/i, ^^b^t, ^^b^^t in the parameter spaces of ei and m-si for different values. 

common region of these two quantities. Existence of solutions satisfying both Ar^re^ and 
Arxefi signifies success of 3+2 model to explain the observed lepton universality violation in 
TT and K decays. Now we turn to Dg decays, whose rates have been experimentally measured 
even though not sufficiently accurate yet, to see if we are able to determine Aro^fir- 

Arxefi and Aro^nr are presented in Fig. M where one notices that within the parameter 
ranges: ei G (0, 6 x 10"^), £2 G (0, 5 x 10^^), m,i G (1,140) MeV, and m,2 G (0,500) MeV, 
both ArKe^J. and Avd^^t- can be satisfied. Especially, in these parameter ranges Ar^e^t and 
AvKe^i are also satisfied. 

The common region for Ar^e^t and Ar^^f^T is shown in Fig. W\ 

The similar analyses have been carried out for D, B, and Be mesons. Due to shortage 
of data to calculate corresponding Ars, we adopt the parameter ranges obtained by fitting 
Ar^re^, ArKe/j. and Ar^^iiT, to the leptonic decays of D, B and Be and investigate if there 
exist common region for Arp^.^ — Arp^j. (P = D^ B, or B^). The results are shown in Fig. 
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FIG. 6: The common solutions of Ar^^e^ and Ar/5^^^ in 3+2 scenario. 
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FIG. 7: The common region of Ar^g^ and Ar/j^^,- in 3+2 scenario. 
[8] and Fig. [9} 

V. TESTING LFU AT BESIII 

Experimentally, measurements on pure leptonic decays of D and Dg mesons have been 
carried out by many collaborations: via e^e~ annihilation at Z^ mass pole [261428] . at 
T(45) mass [291 [30] , and at v^ = 3.773, 4.040 or 4.170 GeV [SH [32], respectively. To test 
violation of the lepton flavor universality, very high accuracy is necessary. But most of the 
the aforementioned experiments suffered from high background contaminations, so do not 
meet the high accuracy demand. In this aspect, the electron-positron colliders prevail over 
others. Because charmed mesons are produced in pairs, one can accurately measure the pure 
leptonic decays based on the double tag method. 

For example, the e^e^ annihilation experiment around the 3.773 GeV, where just above 
the DD production threshold, a charmed meson and its anti-particle are produced in pair, 
i.e. iIj{3770) — ?■ DD. If one fully identifies D in one event, which is called as a singly tagged 
D meson, there must exist a D meson in the recoiling side against the tagged D meson. And 
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if one reconstructs the whole DD pair in the analysis, the event will be called as a double 
tag event. Thus, in an event which consists of a singly tagged D~ , the pure-leptonic modes 
can be selected from the final states of D^ decays, and the absolute branching fractions 
would be well determined. 

For the measurements around 4.040 or 4.170 GeV, situations are not much different except 
DD being replaced by D^Dg or D^D*^^ + c.c 

At the BES III experiments, charmed mesons are collected at 3.773 and 4.040 GeV, 
respectively. Here we present a Monte Carlo (MC) simulation at this two energy points 
to discuss the experimental sensitivities of searching for pure-leptonic decays that can be 
reached in the future. 

The MC events are generated with the BES III offline Software System |33j, where the 
particle trajectories are simulated with the GEANT4 [34J based package [35] for the BESIII 
detector [36] at the BEPC-II collider. 

The events used in this discussion are generated as e"^e~ — )■ ^(3770) — )■ DD and e"^e~ — >• 
^'(4040) — )■ DgD^ at the cm. energy y/s = 3.773 and 4.040 GeV, respectively, where the 
DD and D^Dg mesons are set to decay into all possible final states with the branching 
fractions cited by PDG[22]. 

Totally ~ 1.23 X 10^ DD and ~ 6.20 x 10*^ D^Dg events are generated at y/s = 3.773 
and 4.040 GeV, respectively, corresponding to an integrated luminosity of ~ 20 fb^^ ■?/'(3770) 
and 'ijj{AOAO) data assuming a'^^ = 6.14 nb[2Z] and cr°^+r,- = 0-31 nb[3H], which contains 
~ 7.2 X 10^ D^D^ pairs, ~ 5.1 x lO'^ D+D' pairs and 6.20 x 10^ D^D^ pairs respectively. 
The BEPC-II collider is designed to work with an instantaneous luminosity of 10^^ cm~^s~^ 
around '?/^(3770). As a conservative estimate, a data sample with an integrated luminosity 
of about 20 fb^^ can be collected during more than 10 years' running. 

The singly tagged D~ and Dg events are reconstructed in 9 hadronic decays of 
D- -^ K+7i-7v-{50%), D- -^ K'^sn-{52%), D' -^ K^gK-{A8%), D' -^ K+K-n-{40%), 
D- -^ K+TT-TT-n%28%), D- -^ tt+ n- n- {5Q%) , D' -^ K^-n- -n^ [27%) , D' -^ 
K+Ti+Ti-Ti-Ti-{21%), D- -^ is:^7r-7r-7r+(31%) and 9 hadronic decays of D^ -^ 
/s:^fs:-(46%), D^ -^ K+K-7i-{39%), D^ -^ K+K-n-n^ (12%), D^ -^ /s:^ir+7r-7r-(24%), 
Dg -^ 7r"7r+7r"(52%), Dg -^ n'r], eta -^ 77(41%), Dg -^ n'r]', Vj -^ TX+Ti^r], r] -)■ 77(21%), 
Dg — > -K^Vj ,v( — > 7p''(34%), D^g — )■ p"?7(17%) constituting approximately 29% of all D~ 
decays and 30% D~g decays, respectively, where the numbers in brackets are reconstruction 
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efficiencies. 

Tagged D^ and D~^ events are selected by two kinematic variables based on the principles 
of energy and momentum conservations: (1) Difference in energy 

AE = Ei-Ek, (20) 

where Ef is the total energy of the daughter particle from D~ or Dg in one event and i?b 
is the e~^/e~ beam energy for the experiment, is recorded to describe the deviation from 
energy conservation caused by experimental errors. (2) Beam-constrained mass 



Mbc = V^b - i^^t^) (21) 

is calculated to reduce an uncertainty brought by experimental errors when measuring the 
momenta of the produced particles. By this definition, the energy Ef in the expression of 

Ml,,^Ef-p', 

for the D invariant mass is replaced by E\^ = i5c.m./2, where Ec.m. is the cm. energy at 
which D^D^ pair is produced. 

The total energy and momentum of all the daughter particles in D^ and Ds decays 
must satisfy the Energy Conservation (EC) principle, generally one needs to introduce a 
kinematic fit, including energy and momentum constraints and some relevant corrections, 
to reject those not satisfying EC, but being recorded due to an uncertainty of experimental 
measurement. This replacement of the real invariant mass by Mbc partly plays the role. 

Moreover, events are rejected if they fail to satisfy the selection constraint \AE\ < 3 x o"ae, 
which is tailored for each individual decay mode, and (Jab is the standard deviation of the 
AE distribution. 

As the D^ and Ds events are correctly tagged, a peak in Mbc spectrum would emerge 
at the position of D~ or D^ mass. Finally, if there are more than one combinations in one 
tagged event, the one with the smallest \AE\ is retained. After considering the detection 
efficiencies of each tag mode, 10837045 ±6122 and 549811 ± 1593 tagged D~ and Ds events 
are obtained based on data samples of about 20 fb^^, respectively. 

At the recoiling side against the tagged meson, the other charmed meson decays into 
a charged lepton and a neutrino. Since the neutrino does not electro-magnetically interact 
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with detector matter, it cannot be recorded in the detector and contributes a missing energy, 
therefore, a kinematic quantity is defined 

MLs = (^b - Ei^f - {-to-/D-, - ti^)^ (22) 

where P d- idi i^ ^^^ three-momentum of the fully reconstructed D^ /Dg, and Ei+ {pi+) 
is the energy (momentum) of the candidate lepton. The spectrum of M^^^g for the signal 
events should produce a peak near zero because the neutrino mass is very tiny. 

To select the purely leptonic decay event of D^ — )• I'^ui, only one charged track identified 
as electron/muon/pion and no isolated photons are allowed at the recoiling side. 

For the rare process of D~^ — )■ e~^i>e, the signal number is determined to be 1 by counting 
the signal window of M^jgg. We set an upper limit on the number of signal events for 
D^ — 7- e+z/g to be 4.36 by using the Feldman- Cousins method [39] in absence of background 
at 90% confidence level. The upper limit on the branching fraction for D^ — )■ e'^Ug is 
B{D+ -^ e+Ue) < 8.5 x 10"^ 

For the decay of D^ — )• /i^z/^, the number of simulated signal events is obtained to be 
^D+^a+u, ~ 2611.1 ± 55.0 by fitting the M^j^g distribution, and the efficiency for recon- 
structing D^ —7- yU^z/^ against the tag side is estimated to be iV^*^ +^ = (63.82 ± 0.15)%. 
Therefore, the branching fraction for D^ — t- /i"*"//^ is calculated to be B{D^ — )■ /U+z/^) = 
(3.74±0.08(stat.)) x 10"^ 

For D^ — )• T^u^, with r+ — )• vr+I/^, the missing mass squared M^^^g for the candidate 
events is calculated as 

M^jgg = (Eb - E^^f - i-t^- - t^^y. (23) 

Since there are two missing neutrinos, the M^jgg does not have a narrow peak as in the 
decay of D~^ — )■ fi^u^. To suppress the background from D^ — )• Kln'^ by missing K^ 
and D^ — t- fi^i'^ by /i/tt mis-identification, the M^^gg is studied in two cases defined by 
^EMC ^ Q 2 g^j^j-^ ^EMC ^ g_2 QgV, resulting the number of signal events for D^ — t- t^Uj. to be 
312.1 ±28.2 and 242.9 ±26.1, respectively. With these signal events, inputting the detection 
efficiency of (17. 16 ±0.14)% and (13.33 ±0.12)%, the branching fraction is determined to be 
B{D+ -^ T+Vr) = (1.68±0.12(stat.)) x 10-^ which is averaged by (1.68±0.16(stat.)) x 10^^ 
and (1.68 ± 0.18(stat.)) x 10^^ for the two cases. 

Following the similar analysis, the number of signal events for D^ — t- Z+z/; is determined 
to be N"Jf + = 1, N^Jf ^ = 2335.5 ± 55.1 and N""^^ ^ = 11484.1 ± 110.2, with 
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the corresponding detection efficiencies to be e°^+ ^ = (46.94 ± 0.16)%, e^+ ^ = 
(70.41 ± 0.14)% and e°^_f _,_ = (39.41 ± 0.16)%, respectively. Inserting the numbers of 
events and upper hmit at 90% CL. for Dg — > e"'"z/e, the branching fractions are calculated 
to be B{D^ -^ e+z/e) < 1.69 x 10-^ B{D^ -^ /i+z/^) = (6.03 ± 0.14(stat.)) x 10-^ and 
B{D^ -^ T+Vr) = (5.30 ± 0.05(stat.))%. 

Based on these simulation results, the ratio of decay rates to different leptons can be 
obtained to test the lepton universality. With 20 fb~^ BES III data samples at cm. energy 
of 3.773 and 4.040 GeV, the ratios are 



and 



^^^- = T(D+^T + u) ^ °-^^^ ^ 0.017(stat.) (25) 

for charm meson, and 

RD,,r = r' n! 3* ^ + '''■' = 0,114 ± 0,003(stat,) (27) 

for the strange-charmed meson. The experimental sensitivities for the above measurement 
with 10 years' (20 fb"-*^) data accumulation are A r/p ~ 7.35% and A t/m ~ 2.50%, however, 
the size of these huge data samples cannot present a solid estimate for the electron decays. 
The theoretical expectation for electron mode is at 10^^ level, to challenge this limit, there 
should be a desperate running time for BES III experiment. Fortunately, it will not be a 
problem if one can build a r-charm factory with an increasing of the luminosity of about 
100 times. 

The luminosity of BEPC II is not high enough, so that to fulfill the job, one needs at 
least 10 years machine running. However as suggested, the planned charm-tau factory will 
greatly enhance the luminosity and with the new facility, we expect that in a few months a 
sufficiently large database could be collected and then one may have required accuracy to 
testify the lepton universality. 
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VI. DISCUSSIONS AND CONCLUSIONS 

In this work we study the lepton universahty in the 3+1 and 3+2 scenarios. The analy- 
ses indicate that by adding one sterile neutrino to the SM, i.e., the 3+1 scenario, the Ars 
reflecting the differences between experimental data and SM predictions in the leptonic de- 
cays of various pseudoscalar mesons cannot be accommodated simultaneously for various 
mesons. Therefore, the 3+1 scenario is attractive for its simplicity, but does not meet the 
data. Whereas the 3+2 scenario has less tension in accordance with the data. However, 
although there still exists difficulty in choosing proper a mass range for the second sterile 
neutrino, the current experimental data on Ar^g^, Ar^e^, and Ar£,^^j. can be well accom- 
modated. This result motivates one to incline to involving more sterile neutrinos, i.e., the 
3+3 model. With the parameter ranges selected by fitting the data for vr, K and Dg, the 
predicted values of Ars for D, B, and Be are obtained which will be tested by the future 
experiment such as the Z-factory and LHCb. 

Usually it is believed that the ft — r symmetry [ID] holds at high-energy scales, but breaks 
during the evolution to low-energy, so that the 3x3 PMNS mixing matrix for active neutrinos 
deviates from the original symmetric textures. The violation of lepton universality, especially 
that for fi — T universality, might be a low-energy behavior as the universality precisely holds 
at high-energy scale, say, the GUT or the see-saw energy scales. With these assumptions, 
it motivates one to relate the breaking of the fi — t symmetry to the violation of the fi — t 
universality. The idea is that the symmetry breaking leading to real PMNS matrix and the 
LFU violation may originate from same source and caused by the smae mechanism. Thus 
both of them serve as the low-energy manifestations of the symmetry beaking. 

In order to obtain a negative value for Ar, Rp^a should be smaller than the value predicted 
by the standard model R^p- Then one can obtain 

'Zfjf \Uo,,i+3\'^iml + ml-){ml -ml- mli)X{ml,ml,m'^J ' 

i.e., the ratio of the contributions from sterile neutrino should be larger than the SM pre- 
diction i?SM^. 

Obviously, ignoring high order QED radiative corrections, only the mass of the concerned 
pseudoscalar meson enters the game, but not the identities of its constituents, thus we can 
relate Ars of various mesons to each others. This conclusion is viable for checking the 
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scenarios discussed in the introduction. 

Our numerical analyses indicate that checking violation of LFU presents a rigorous re- 
quirement for our colliders, therefore a high-luminosity charm-tau factory and/or B-factory 
are necessary to draw solid conclusion. 
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